Analytical Solution of the Proprio-Graviceptive equation for shoot gravitropism of 

plants. 



Renaud Bastienj^'^'-^^^'j^] Bruno Moulia,^''^ Stephane Douady,^ and Tomas Bohr^ 

^Institut Jean-Pierre Bourgin, UMR1S18 INRA-AgroParisTech, 78026 Versailles, France 
"^INRA, UMR 547 PIAF, F-63100 Clermont Fd Cedex 01 
^Clermont Umversite, Universite Blaise Pascal, UMR 547 PIAF, BP 10448, F-63000 Clermont-Ferrand 
'^Matiere et Systemes Complexes, Universite Paris- Diderot, 
10 rue Alice Domont et Leonie Duquet, 75025 Paris Cedex 13, France 
^Center of Fluid Dynamics and Department of Physics, 
Technical University of Denmark, 2800 Kgs. Lyngby, Denmark 

We derive the analytical solutions to the second order generalised gravi-proprioceptive equation 
given in our recent paper [I]. These equations show how plants adjust to the surrounding gravitation 
field and highlight the fact that the plant must be able to not only sense its local posture with 
respect to the gravitational field, but also to sense its own local curvature. In [l] we obtained 
explicit analytical solutions of these equations in terms of (sums of) Bessel functions, and in the 
present paper we derive these solutions. 



Gravitropism is a slow movement by which plants re- 
orient their growth in response to gravity. In a recent 
paper [T], we have studied the gravitropic response of a 
broad range of plants by placing them horizontally in a 
dark room and monitoring how they curve upwards in 
response to gravity. It was shown that the minimal dy- 
namical model accounting for these observations is not 
the usual gravisensing law, where the plant simply senses 
the local tilt angle, but what we have called a gravi- 
proprioceptive law, where one takes into account that the 
plant can, in addition, sense its own local curvature. In 
that paper, exact solutions for the linearized versions of 
both the graviceptive and the gravi-proprioceptive law 
were given and compared to experiment. We here give 
the derivation of these solutions 

The plant is described as a slender rod with arc length 
s going from s = at the base to s = L at the tip. 
The local tilt angle A{s,t) describes the orientation with 
respect to the vertical gravitational field, such that yl = 
is vertical. The generalized (gravi-proprioceptive) law of 
gravitropism describes how the local curvature C(s,i)) 
changes in response to the local tilt angle and the local 
curvature itself as 



dC{s,t) 
dt 



-PA{s,t)~^C{s,t), 



(1) 



where the curvature and the tilt angle are related by 



C{s,t) = 



dA{s,t) 
ds 



(2) 



and where the coefficient j3 controls the graviceptive term 
(sensing of tilt angle) and the coefficient 7 controls the 
gravi-proprioceptive term (sensing the curvature). 

The subject of the present paper is to solve this equa- 
tion with the experimentally relevant boundary and ini- 
tial condition. Thus we assume that the plant is initially 
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straight with an angle, say Ai from the vertical. Since 
([1]) is linear in A, we can, without loss of generality, take 
Ai — 1. The solution with a different Ai will simply 
be obtained by multiplying our solution by A.^.Thus we 
assume the initial conditions 



A(s,t = 0) = 1 



and 



C(s,i = 0) = 



dA{s,t^ 0) 



^ 



(3) 



(4) 



for all < s < L. Further we assume that the plant is 
clamped in the sense that tilt angle at the base remains 
equal to Ai = 1 for all future times. Thus the boundary 
condition is 



A{s = 0,t) = 1 



(5) 



for alH > 



I. THE GRAVICEPTIVE MODEL 

We first consider the purely graviceptive model, i.e., 
^ with 7 = 0: 



(6) 



To solve this equation, we first set P — 1 with out loss 
of generality simply by using the scaled time pt. Now, 
using that 

dA{s,t) 



C{s,t) 



ds 



(7) 



we can write ([6]) as 



d^A{s,t) _ 

dsdt ^ 

where we use the short hand notation d/dx = dx- Note 
that C satisfies the same equation. 
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A. Solution by separation 

We shall now try to find separated solutions of 
From our numerical work it seems that the solutions de- 
pend on y/s and ^/i and since both s and t are assumed 
positive, we can introduce the new independent variables 
(C,?7) as 



s = — and t = S^rj 



1/2 



(9) 



(10) 



and obtain 



ds = m) + (5,77) 5, = ^vd^ ~ ^v^C' d„ (11) 
dt = (m + (dtV) dr, = Iv-' + dr, (12) 



Using this, we find 
1 



4^2 



2^)2 N 

vn/ 



(13) 



Note that the cancel and that the invariance of (|6| on 
exchanging s and t implies an invariance when 77 — >■ r]^^. 
In these variables, ([S]) for an arbitrary function ^1/ can be 
written 

+ - vdv + V^d^,, + 4e') V) = (14) 

If we now assume that the solution can be written in the 
separated form 



*(c,'7) = /(^)5(e) 

we get the two eigenvalue equations 

cV(e) + e5'(0 + 4evo-Ag(e) 



Substituting y — logry gives (17) the form 



with solutions 



with 



fiy) 



2p + X = 



(15) 



(16) 
(17) 



(18) 
(19) 
(20) 
(21) 



and then ( 16 ) becomes 



e9"i0 + e g'iO + i^e - P{2 - p)) giO = (22) 



If we finally substitute a; = 2^, this becomes 

x^g"{x) + X g'[x) + {x^ ~ p{2 - p)) g{x) - (23) 

which is recognised as Bessel's equation of order n, where 
m? = p{2 — p). For n to be a real number we must restrict 
p to lie in the interval p G [0, 2]. In this interval we can 
write the full solution as 



p/2 



(24) 



If n is to be an integer we must choose p = 0, 1 or 2 
giving 71 = 0, 1 and 0, respectively. To find a solution for 
A, we must satisfy the boundary condition A{s — > 0, t) = 
A{s, t 0) = 1 and since 



(25) 



for small st, we must choose n = p = 0, i.e., 

A(^,77) = Jo(20 (26) 
or, returning to the variables s and t and re-inserting /3, 

A(s,t) = Jo(2v/^) (27) 
We can also find the curvature 



C(s,t) = dsA{s, i) = y y Ji (2v/^) 



(28) 



which, since C also satisfies (|8|, also has the form (24 1. 
For small s, the s-dependence cancels since Ji{x) Rix/2 
for small x, and we get 

C(s,i)-^/3t (29) 

for small s or t and thus C(s, i — > 0) = as we demanded. 

B. Solution by Laplace transformation 

We could have solved ([s]) more directly Laplace trans- 
forming is (in s), which gives 



where 



pdtA - dtA{s = 0,t) = -M, 



A{p,t) = / A{s,t)e-P'ds 
Jo 



(30) 



(31) 



is the Laplace transform of A{s,t). From the boundary 
condition ([5| we get dtA{s = 0,t) = and 



dtA^~'-A, 



(32) 



which has to be solved with the initial condition A(s, t 
0) = 1 ^ A{p, t = 0) = 1/p with the solution 

i 1 _ 3t 
A = —e p 
P 
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We see that A{p, t) is analytic in the complex p-plane 
except at the isolated (essential) singularity p = 0. The 
inverse Laplace transform A(s, t) can simply be found as 
the residue (see e.g. |^ Theorem 8.2.1) 



A{s,t) 



2m 



-dp = Jo(2V/^) (33) 



where the last equality is a special case of the well-known 
identity derived e.g. in [3]: 



1 



27r« Jn{o,e) P' 
valid for any n > — 1 



' P / s \ "/■^ 



(34) 



with the solution 
A ' 



1 i+Pt 



l+p p{l+p) 
1 /I 1 



I + P \p l+p 



e e p 



I _.t 1 

-e e p H 

p l+p 



l-e-'e^p), (43) 



where A(p, t = 0) = l/pin accordance with the boundary 
condition. It correctly approaches the stationary state 
since 



i(p,t) 



1 



l+p 



An 



(44) 



II. THE PROPRIO-GRAVICEPTIVE MODEL 

We now consider the generalised equation of gravit- 
ropism 

dtC{s,t) = -f3A{s,t)--fC{s,t) (35) 

and with scaled time T — "ft and S ~ {j3 /"f)s we obtain 

dl^A = -A- dsA (36) 

In the following we drop the capitals and our equation 
becomes 



di^A =-A- dsA 



(37) 

again with the initial conditions ([3| and Q and the 
clamped boundary condition ([s]). 

Note that (37) has that stationary (asymptotic) solu- 
tion 



Ao{s) 



(38) 



satisfying the boundary conditions [3] and |4] One might 
try the substitution yl(s,i) = B{s,t)e~'^, since the equa- 
tion for B would simply be Q. However, the boundary 
condition ([sj would be B{s = 0,t) = e* which is time 
dependent and complicates matters. 
The Laplace transform is 



A{p,t) = / A{s,t)e-P'ds 
Jq 



(39) 



and the Laplace transform of 37 (with "/ — {3 — I) gives 

pdtA - dtA{s = 0, t) = -i - {pA ~ A{s = 0, t)) (40) 

From the boundary condition, we know that A{s — 
0, t) is constant in time and thus dtA{s = 0,t) = 0. Also 
A{s = 0,t) = I so we obtain 



pdtA 



-A-pA + l 



(41) 



for t — )■ oo when p > and A(p, t) is analytic in the 
complex j3-plane except at the isolated singularities p = 
and p = — 1. 

Thus the inverse Laplace transform A{s,t) is simply 
the sum of the residues [2] of the function 



A{s,t) 



E 



1 



2711 Jn(„, 

singularities p* 



e e p 



l+p 



1 — r — — 

1 — e e p 



(45) 



at each of its singularities in the complex p-plane. The 
pole at p = —1 is a simple pole, but the residue is zero 

— i — - 

since the numerator 1 — e~ e p is zero for p = —1. We 
therefore only need the residue at p = 0, which is an es- 
sential singularity. Thus we have to evaluate the contour 
integral of e^*v4(p, i) over a closed curve S1(0, e) encir- 
cling the origin. The contribution of the term e^'' /{I +p) 
vanishes, since it has no pole at p = 0, and we can write 



where 



and 



A{s,t)=Ai{s,t) + A2is,t), 



g-t y gsp-p 

Alis,t) ^ ^ f dp 

27r« 7o(o,e) P 



27ri /o(o,<r) l+p 



where first term is given by p3[ ) 

Ai(s,t) = e"*Jo(2Vst). 



(46) 



(47) 



(48) 



(49) 



Near the origin (more precisely, as long as | p |< 1) we 
can expand 



dtA+ -A ^ - 

P P 



(42) 



1 



l+p 



l-p + p^ 



(50) 
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so that 



write 



A2(s,i) 



2m 



E 

n=0 



e'P-vi-pYdp 



0.(0,, 



(51) 



n=0 



where 



2TTi 



n(0,e) 



e''P-p(_p)«dp. (52) 



We define z = — l/p and thus dp = dz/z^, whereby 



2m 

t 

s 



Sl(0,e 
+ 1 



J„+i(2Vsi) 



(53) 



for any n > Q. Note the sign change since the contour f2 
is now traversed in the clockwise direction. Thus we can 



m— 



E 



J„(2VsO 



(54) 



and the full solution is 

^/2 



m=0 



^(-1)"(-) J„(2v^) (55) 



which, going back to the original variables using t ^ jt 
and s — >■ {P/j)s gives 



2+\ 2 



A(s,t)=e-^*5] ( ^) J„(2V/3sO 



/3s 



m=0 

= e-^^/^-^(-l)" 4- J„(2y^) (56) 
,1=1 

If the initial angle is Ai, we get the same solution, but 
with the factor Ai multiplying the solution ( 56 ) . 
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